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A new theoretical method is introduced to study coherent electron transport in an interacting 
multilevel quantum dot. The method yields the correct behavior both in the limit of weak and 
strong coupling to the leads, giving a unified description of Coulomb blockade and the Kondo effect. 
Results for the density of states and the temperature dependent conductance for a two-level dot are 
presented. The relevance of thessrfesults in connection to recent experiments on the Kondo effect 
in semiconducting quantum dotsocl is discussed. 
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The proper description of coherent electron trans- 
port in the presence of strong electron-electron inter- 
actions has been one -of the central issues in the field 
of mesoscopic systems!]. Semiconducting quantum dots 
(QDs) provide an almost ideal system where the pre- 
dictions of the theory can be tested and new effects 
could be searched for. For instance, recent experi- 
ments have demonstrated the possibility of exploring the 
Kondo effect, a prototypical correlation effect, using this 
technologyBB. 

From the theoretical point of view, correlation effects 
have been mainly analyzed in this kind of systems by 
means of the simple Anderson model with a single spin- 
degenerate leveloO. Within this model the Kondo ef- 
fect arises due to fluctuations in the spin of an unpaired 
electron!! So far, there has been few attempts to include 
a multilevel spectrum for describing either single QDs 
with quasi-degenerate levels or coupled QDs beyond a 
rate equation approacfjj. The actual situation in semi- 
conducting QDs should require the inclusion of a multi- 
level spectrum whenever many single-particle states with 
a small level separation are involved in the transport pro- 
cess. This could be the case in experiments measuring 
the phase of the transmission amplitude through a QD 
in the Coulomb blockade (CB) regimeS and also in recent 
experiments on the Kondo effectEm. 

The aim of this letter is to introduce a new theoretical 
approach for describing correlation effects in multilevel 
QDs. The approach is constructed to yield the correct 
behavior both in the limit of infinite and vanishing charg- 
ing energy. This is achieved by introducing an interpola- 
tive self-energy for the one-electron Green functions, an 
approach which has been successfully applied to several 
interacting systems including the equilibrium and non- 



equilibrium Anderson modelsSi. This type of approach 
has been recently rediscovered and applied to analyze the 
M ott Jxansition in the Hubbard model for arbitrary band 
fiwngliM 

For describing the multilevel QD we consider a model 
Hamiltonian which is a generalization of the single-level 



Anderson model, H = H, 



dot 



leads 



Ht , where 



H. 



dot 



E, 



uncoupled QD 



U J2i> m n m ni corresponds to the 
= d] n d rn )\ #l cac ls = J2keLM e kc\;Ck 



to the uncoupled leads, and Ht = E 



m,k£L,R 



h.c. describes the coupling between the dot and the leads. 
The labels m and I in H denote the different dot levels 
including spin quantum numbers. 

Our main objective is the determination of the dot 
retarded Green functions from which the different level 
charges and the linear conductance can be obtained. In a 
frequency representation they can be written as G m (u>) — 

[u> - e^ F - E TO (w) - iT rn . L (uj) - iT m . R (uj)] , where 
e^ F = e m + U J2i^ m 711 ^ s tne Hartree-Fock level (we 
adopt the notation ni =< hi >), iF m T(R){^) — 
t 2 m k / (to — gfe + i0 + ) are the hybridization terms 



eL(R) b m,k 



Efc 

describing the coupling to the leads and S m (w) is a self- 
energy that takes into account electron-electron interac- 
tions beyond the Hartree-Fock approximation. We-,shall 
neglect the indirect coupling between dot stateslla and 
adopt the usual approximation of T™^, T„ lt n being inde- 
pendent of energy. In our approachErcl, we look for an in- 
terpolative self-energy yielding the correct U/T m — > 0, oo 
limits. 

In the U/T m — > oo or "atomic" limit, G m (col can be 
obtained using the equation of motion methodic, which 
yields 
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w — (m — U + i0 + 
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(1) 



where M denotes the total number of one-electron lev- 
els in -H^ot' ^ n this expression all possible charge states 
of the dot give a contribution. Their evaluation requires 



the knowledge of many particle correlations < h\h% >, 
< h\hihz >, etc. However, for sufficiently large U fluc- 
tuations in the dot charge with respect to the mean 



1 



charge N by more than one electron become negligible 
and Gm (w) is accurately given by the three poles ex- 
pression 

4' 
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uj - e m - U(N - 1) + i0+ 
4" 



£W + zO+ uj 



U(N + 1) + i0+ 



(2) 



where TV = /nf [A/]. In order to yield the exact first three 
momenta of Eq. (1) the weight factors A 7 ^ should satisfy 
the following sum rules 
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N-l 



A m 4- A m — 1 



(N - + NA™ + (N + 1)A% +1 



(N - 1) 2 ^_! + iV 2 ^ + (N + 1) 2 A™ +1 = 

ni+ < hh > m , (3) 

where < hh > m = J2(i^k)^m < n i n k >■ The special case 
N = (JV = M - 1) has to be treated as N = 1 (N = 
M — 2). Notice that this expression for Gm''(w) is fully 
determined by the average charges ni and the two-body 
correlations < fiihf. >. From Gm (w) one can define an 



a m U 2 (uj ~e m + i0+) + b m U 3 



(uj - e rn + i0+) 2 + c m U{uj - e m + i0+) + d m U 2 ' 



(4) 



where a m — (AT - n m ) [1 - (J\f — n m )] + < hh > m ; c m = 
Af-n m -3N; d m =< hh > m +3N 2 -l-(3N-l)(jV-n m ) 
and b m = N 2 (l — N) — (j\f — n m )d m . On the other hand, 
in the f /Cm ~~ * limit E m is given by the second order 
expressionO 



Sg'H = U 2 J2 de x de 2 dt 



Pm{ei)pi(e 2 )pi(e 3 ) 
uj + e 2 - e\ - e 3 + i0+ 
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"atomic" self-energy, S 

Using Eqs. (2) and (3) it can be shown that Em can be 
written as the ratio of two polynomials in uj of the form (4) for E^ 1 ' , which yields 



where /j = /(ej) denotes the Fermi distribution func- 
tion at the leads, and p m (uj) are effective densities of 
states given by np m (uj) = Y m /((uj - e m ) 2 + Y 2 m ). The 
effective levels e m are introduced in order to fulfill the 
Fermi-liquid relations associated with charge conserva- 
tion (Friedel sum rule (FSR)EJ) as discussed below. 

In order to determine an interpolative scheme between 
the two limits let us first notice that — > U 2 a m / (uj — 

e m ), where a m = J2i^ m "K 1 _ "0> wnen r m/ w °- 0n 
the other hand, a m — > a m in the small U limit and thus 

Em — > U 2 a m /(uj — e m ) in this case. These properties 
suggest to define the interpolative self-energy, replacing 
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e m by a m /S™ + Ae m , where Ae r 



i(2)' 



a m + 



(a m Ae m /f/ + 6 m /a m ) (e^E/) 



1 + (2Ae m /U + c m ) (Y$/a m U) + (c m Ae m /U + d m ) /a m U 



in Eq. 



(6) 



This expression provides the generalization to the mul- 
tilevel case of the interpolative self-energy first intro- 
duced ina for the single level Anderson model. The single 
level case (M = 2) is readily obtained from Eq. (6) when 
N = 1 taking < hh > m = 0. It is evident by construction 
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the small U limit. 

The final step in this approach is to determine the level 
charges n m , the correlation functions < hih^ > and the 
effective levels e m selfconsistently. The charges and the 
correlations functions are determined through the rela- 
tions n m = — f (uj)ImG m (uj)duj /tt and 

< hih m >= — / f(uj)ln\[T, m (uj)G m (uj)]duj. (7) 

Eq. (7) is an exact relation that follows from the equa- 
tion of motion of the retarded Green functions. It is im- 
portant to stress that the self-consistent determination of 
the two-body correlations is essential to get the correct 
values of the charge for large U. The effective level is 
determined by imposing the condition: 



/Mini 



G m (uj) 



<9E m 
duj 



duj = 0, 



(8) 



which at, zero temperature reduces to the Luttingcr 
theoremEj ensuring the fulfillment of the FSRli-J n m — 
— Im[ln G m (EF)]/TT . In RefB we showed that the con- 
dition of consistency between the effective and the final 
charges is nearly equivalent to imposing the fulfillment 
of the FSR for the simple Anderson model. Simila r self - 
consistency conditions have been proposed in Refs.EM. 

We have first applied this formalism to the case of a 
doubly degenerate level (four-fold degeneracy when in- 
cluding spin) which is a simple generalization of the sin- 
gle level Anderson model. The inset in Fig. 1 shows 
the charge per level Tlffi clS cl function of the leads Fermi 
energy Ef- For the case shown in Fig. 1, correspond- 
ing to e m = 0, T m /U = 0.075 (same for all levels) and 
zero temperature, one can observe a modulation in the 
charge resembling the typical Coulomb staircase of the 
^m/U — > limit (shown as dotted line for comparison). 
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FIG. 1. Density of states for a fully degenerate two-level 
dot with Tm/U — 0.075 and e m = for different values of the 
leads Fermi energy Ef- The inset shows the charge per level 
and spin as a function of Ef (full line) as compared to the 
r m /(7 — > case (dotted line). Energies are measured in units 
of U. 

Fig. 1 also shows the density of states (DOS) associ- 
ated with the dot levels for different values of Ep/U — 
0,0.5, 1.0, 1.5 corresponding to the steps and the center 
of the plateaus in the charging curve. For Ep/U = 0, 
n m ~ 0.125 (the total charge on the dot being ~ 0.5e) 
the system is in the so-called mixed-valence regime. The 
spectrum is in this case similar to the one found us- 
ing large iV-expansions in the U — > oo limitE3 show- 
ing a resonance at a renormalized level just above Ep. 
For Ep/U = 0.5 and 1.5 there are approximately one 
and two electrons respectively inside the dot. For these 
cases the system is in the Kondo regime, the DOS ex- 
hibiting a sharp peak around Ep (Kondo resonance) half 
way between two broader resonances (of width ~ T) ap- 
proximately separated by U, corresponding to the ele- 
mentary charge excitation energy. It is worth noticing 
that the weight of these resonances is the same for the 
Ep/U — 1.5 case due to the electron-hole symmetry in 
the half- filled dot, while there is a pronounced asymme- 
try for Ep/U — 0.5. For a half-integer occupation of 
the dot like in the Ep/U — 1.0 case the DOS exhibits 
a more complex structure with three broad resonances 
around E — 0, U and 2U, and a narrow Kondo peak still 
present at Ep. The overall shape is reminiscent of an 
average between the cases with one and two electrons in 
the dot. 

The evolution of the Kondo peak as a function of Ep 
and temperature should be reflected in the dot linear 
conductance, G, which caa be readily obtained from the 
Green functions (see Ref £3) . The conductance as a func- 
tion of Ep and different temperatures is shown in Fig. 
2. When approaching zero temperature the conductance 
behaves like G = 4(e 2 //i) sin 2 [irn m ] as expected from 
the FSR. The conductance decreases very rapidly with 
temperature in the region 0.5 < Ep < 2.5 where the 
DOS at Ep is controlled by the Kondo peak. Outside 
this region one can notice a slight increase of conduc- 



tance with temperature. At temperatures large enough 
to be above the Kondo temperature, T%, which can be 
estimated by the condition G(Tr-) ~ G(0)/2&l, the con- 
ductance tends to exhibit the usual CB peaks at the 
charge degeneracy points. We should remark that our 
approach does not provide an accurate estimate of Tk as 
a function of the model parameters as the exponential de- 
crease of the Kondo peak weight for very large UtB is not 
strictly recovered. This limitation does not affect, how- 
ever, the qualitative behavior of the conductance except 
for a rescaling of the temperature values. 
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FIG. 2. Total conductance for the same case of Fig. 1 as 
a function of Ef for different temperature values T/U : a) 
0.0005, b) 0.0025, c) 0.005, d) 0.01 and e) 0.03. 

So far we have analyzed the case of a completely de- 
generate level. In an actual QD geometrical asymmetries 
could result in an effective splitting AE of the dot levels. 
When AE » T the physical situation could be described 
by a single level Anderson model. It is interesting to 
analyze the evolution from this situation to the quasi- 
degenerate case AE — > previously discussed. This 
evolution is depicted in Fig. 3 and 4 where a splitting 
AE = 0.1U and AE = 0.5U is introduced. Only the first 
half of the charging curve is shown as the total conduc- 
tance is symmetrical with respect to Ep = (3U + AE)/2. 
As can be observed in the charging curves shown as insets 
in these figures, the splitting tends to block the charging 
of the upper levels. This blocking effect is nearly com- 
plete for AE — 0.5U and is already apparent in the case 
of Fig. 3 where AE ~ T. The total conductance at low 
temperatures (case a) is observed to decrease for increas- 
ing splitting approaching the behavior 2(e 2 //i) sin 2 [7m m ] , 
where m denotes the lower level. The overall effect of 
temperature is again to reduce the total conductance in 
the range for Ep where the Kondo effect is present. In the 
case of Fig. 4 the conductance increases with tempera- 
ture in the region Ep ~ 1.5U where the higher level starts 
to be populated. One can also notice for this case that 
the conductance peaks associated with the lower level are 
not completely symmetrical due to the influence of the 
second level. Both effects, the increase of conductance 
with temperature between each pair of peaks as well as 
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the asymmetry, are present in the experimental results 
of Ref.El and can be considered as a manifestation of the 
multilevel structure of an actual QD. 
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FIG. 3. Total conductance for a two-level dot with splitting 
AE/U = 0.1 as a function of Ep for different temperature val- 
ues T/U: a) 0.0005, b) 0.0025, c) 0.005 and d) 0.03. The inset 
shows the charge per spin on the two levels as a function of 
Ep. Only the first half of the curve is shown as the total 
conductance is symmetrical with respect to (3!7 + A)/2. 
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FIG. 4. Same as Fig. 3 for AE/U = 0.5. 

The most remarkable consequences of this multilevel 
structure would be observed in a situation correspond- 
ing to a smaller splitting where "collective" Kondo like 
features like the ones depicted in Figs. 2 and 3 should 
appear. Vertical dots with cylindrical symmetry, as those 
studied in RefJ13, constitute an almost ideal realization of 
the two-fold degenerate case. One word of caution should 
be said, however, regarding the effects of Hund's rule, 
not included in the present approach and which might 
introduce some deviation with respect to the behavior 
depicted in Fig. 2. The results presented in this letter 
will fully apply to a situation where the exchange interac- 
tion between the degenerate levels is much smaller than 



the Coulomb interaction U. 
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